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EXPERIMENTAL EVALUATION OF A PUMP TEST FACILITY WITH 
CONTROLLED PERTURBATIONS OF INLET FLOW 
by William Stevans and Robert J. Blade 
Lewis Research Center 

SUMMARY 

A description of a pump dynamics test facility and the results of its experimental 
evaluation are presented. The facility is designed to experimentally investigate the re- 
lations governing the operation of a pump under periodic fluctuations of inlet pressure 
and flow. A specially designed throttle produces perturbations to the pump inlet flow and 
pressure over a frequency range of 2 to 60 hertz. The perturbed flow entering the pump 
is determined by using the acoustical wave equation to determine the flow in the long in- 
let line and by accounting for the equivalent vibrational flow due to the longitudinal mo- 
tion of the inlet line. The perturbed flow leaving the pump is determined by using a 
specially designed multihole orifice plate. 

Independent checks for establishing the validity of the perturbed data are described 
and applied to the test data. The results indicate that the measurement methods used 
give reliable values for the perturbed flows and that the motion of the test facility must 
be accounted for when calculating the perturbed flow entering the pump. Pipe motions 
of only 7. 62x10 meter (0. 003 in. ) can result in an error of 100 percent in the flow en- 
tering the pump. The difference between the flows entering and leaving the pump was 
identified as the compliance of the pump structure and the water it contained. The larg- 

c o 

est volume difference detected was approximately 0. 06x10” 3 cubic meter (0. 04 in. J ). 

An adaptation of the Frahm Vibration Absorber was used to check the vibration flow 
measurement. The device was found to be effective in reducing the inlet line longitudinal 
motion for each resonant frequency for which it was set. However, it was effective over 
only a narrow frequency band about each resonance point and therefore impractical for 
continuous use over the frequency spectrum covered in the pump tests. 


INTRODUCTION 


A serious problem in rocket launch vehicle performance is the longitudinal struc- 
tural oscillation of the vehicle in its fundamental mode, often referred to as "pogo. ” 
These longitudinal oscillations are a result of the unstable closed-loop coupling of the 
vehicle structure and the vehicle propulsion system. The frequency and time of occur- 
rence of the maximum oscillations are different for each type of launch vehicle. For 
vehicles presently in use, the frequency of these oscillations falls in the range of 5 to 
25 hertz. These oscillations usually occur near the end of the boost phase of flight. 

The unstable loop coupling of the vehicle structure and propulsion system is self- 
sustaining and has the following pattern: The structxiral longitudinal oscillations cause 
variations in the pressures and flows in the propellant feed systems to the pumps; this 
results in a variation in the pressures and flows of the propellants delivered to the en- 
gine combustion chamber and, hence, in a variation in engine thrust. This variation in 
engine thrust produces a longitudinal structural oscillation which completes the unstable 
loop. 

Studies of pogo oscillations indicate that it is difficult to write the accurate mathe- 
matical transfer function for the performance characteristics of the turbopumps if 
there is a variation in their inlet pressures and flows. The analytical study of refer- 
ence 1 represented the turbopump by equations that could be obtained from its steady- 
state operating characteristics. Two recent experimental investigations to define the 
turbopump dynamic operating characteristics are reported in references 2 and 3. In 
these studies certain assumptions or restrictions were placed on the terms of the pump 
transfer function which prevented obtaining an expression in its most general form. 

In order to obtain the pump transfer function in its general form, a experimental 
study of pump dynamics was undertaken at NASA Lewis Research Center. To conduct 
these studies a research test facility was constructed which would allow detailed meas- 
urements of the fluctuations of pressure and flow at the inlet and exit of the pump. 

This report describes the general considerations applied to the design of the test 
facility, gives specific details of its major components, and presents data which demon- 
strate that the necessary measurements to permit analysis of pump dynamic character- 
istics can be made in this facility. 


APPROACH 

In this experimental evaluation of pump performance with perturbations to the pump 
inlet conditions, there are two areas of primary interest. The first is the determina- 
tion of pump resistance and gain terms which are obtained under both noncavitating and 
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cavitating conditions of pump operation. The second is the determination of the compli- 
ance characteristics associated with the pump operating in a cavitating mode. 

For facility design purposes the analytical representation of the pump was con- 
sidered to have the following form: 


P d = (1 + A)P S - BQ d (1) 

This expression relates the perturbed inlet pressure P„ and the perturbed discharge 
pressure P d to the perturbed outlet flow rate Q d> The perturbed pressures and 
flows are periodic sinusoidal variations about the mean value of the respective param- 
eters. The terms A and B are considered to be representative of the pump gain and 
impedance, respectively, and can be complex values. 

When the analytical representation of the pump is of the form of equation (1), there 
are several different methods of establishing values for the terms A and B. A num- 
ber of methods were evaluated and the following is the method used in the experimental 
investigation conducted in this test facility. It was chosen because it permits a solution 
for unique values of terms A and B. Inspection of equation (1) indicates that for non- 
cavitating operating conditions of the pump, where term A is assumed to be zero, 
unique solutions for term B can be obtained by using a single set of dynamic test data. 

For cavitating operating conditions of the pump, both A and B are unknown; this 
precludes obtaining a unique solution for them with a single set of dynamic data. For 
this situation, two sets of dynamic data must be obtained while the pump is operating at 
the same steady- state operating conditions. Rewriting equation (1) in the following form 
indicates a method by which this may be accomplished: 

P d = (1 + A)P g - BQ d = ZQ d (2) 

where Z is the hydraulic impedance at the outlet of the test pump. By changing the 
value of Z, two sets of data in which A and B are unknown can be obtained. This 
allows for the solution of equation (2) by the method of simultaneous equations and the 
obtaining of unique values for A and B. Changes in the value of Z, with the pump 
operating at the same steady-state conditions for both Z values, can be obtained by 
using orifice plates of different pressure drops at the pump outlet. To maintain the 
same steady- state operating conditions of the pump with different orifice plate pressure 
drops, the overall pump outlet pressure drop must be adjusted. This can be accom- 
plished by adjusting a valve, which is placed downstream of the orifice plate, so the 
same overall steady- state outlet pressure drop exists for both orifice conditions. 

The second objective of this experimental evaluation of pump performance with per- 
turbed inlet conditions is the determination of the compliance characteristics associated 
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with the pump when it is operating with cavitation at its inlet. The determination of the 
compliance characteristics is based on the measurement of a compliance flow which is 
the difference between the perturbed flows entering and leaving the pump. This relation 
can be expressed by an equation of the following form: 

% - % = EP s < 3 > 

where E is a measure of the compliance. 

The methods that were used to measure the perturbed flows required in this inves- 
tigation were suggested by experience gained in the line dynamics work at this labora- 
tory (refs. 4 and 5). The use of these measuring methods has a significant effect on the 
configuration of the test facility. 


TEST FACILITY 

The Lewis Pump Perturbations Test Facility, shown schematically in figure 1, is 
a closed-loop system containing approximately 9. 5 cubic meters (2500 gal) of conditioned 
water. The facility is an adaptation to the major components of the Lewis Water Tunnel. 
The components of the Lewis Water Tunnel are the drive system, the water conditioning 
system, and the mean flow measurement system, all of which are described in refer- 
ence 6. The components added for the Perturbations Test Facility are the inlet and out- 
let dynamic flow measurement systems, the inlet and outlet isolation tanks, the inlet 
flow and pressure perturbing system, the inlet line vibration absorber, and the test 
pump. Included as an integral part of these items are the dynamic and steady-state in- 
strumentation. 


Test Pump 

The pump impeller selected for these tests was the unshrouded centrifugal impeller 
which is shown in figure 2. The pump was designed to operate in the Lewis Water Tun- 
nel, and data on its steady-state performance are given in reference 7. 

Dynamic Flow Measurement Systems 

The first flow of interest is the flow entering the pump. This flow is determined by 
measuring the flow in the inlet line approaching the pump and subtracting the equivalent 
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flow due to the motion of the inlet line. The perturbed flow approaching the test pump 
was measured by applying the acoustical wave equation to measurements of pressure 
fluctuations in a long hydraulic line. Detailed derivations of the equations used are 
given in the appendix of reference 4. The use of this method of measurement requires 
that the inlet system be of sufficient length (1) to allow calculation of the flow at the pump 
at low test frequencies, and (2) to allow accurate determination of the acoustic velocity 
of the fluid at the maximum test frequency (see appendix B for a derivation of the method 
of determining the acoustic velocity). The inlet system must also be capable of passing 
the mean flow required by the test pump with a small pressure loss. Based on these 
requirements the inlet flow measurement system is a 18. 3-meter (60-ft) long pipe of 
0. 20-meter (8-in. ) diameter with a 3. 76xl0~^-meter (0. 148-in. ) wall thickness. The 
pipeline is made in four sections which are flanged and bolted together. The flanged 
joints are constructed so a smooth continuous internal bore is maintained in the pipeline. 

The pressure perturbations in the line are measured with five pressure sensors 
equally spaced along the bottom of the pipeline (see fig. 1) with the sensor diaphragms 
flush with the internal bore of the pipe. This placement was made to minimize the 
chance of error in pressure measurement due to gas being trapped near the region of 
the sensor's diaphragm. The number of sensors used allows four independent calcula- 
tions of flow to be made at station 5. 

The second flow of interest is the perturbed flow leaving the test pump. This was 
determined by measuring the perturbed pressure drop across a specially designed mul- 
tihole orifice plate. Data obtained in the line dynamics work (see table I, ref. 5) indi- 
cate that such an orifice plate could be used to measure perturbed flow. 

The design of the orifice plate is based on its dynamic impedance having resistive 
and inertive reactive components. The resistive component is assumed to be the slope 
of the tangent to the steady-state-pressure-drop-against-flow-rate curve for the orifice 
plate holes. The reactive component is a function of the inertial mass of the fluid in the 
orifice plate holes and the equivalent of two Rayleigh end corrections for the holes. The 
effect of the reactance term is reduced to a practical minimum by resorting to a small 
hole diameter (1. 70xl0~^ m, or 0. 067 in. ) and a thin plate (1. 27 x 10 ”^ m, 0 r 0. 050 in. ). 
The difference between the reactance achieved and the theoretical minimum was equal to 
the resolution power of the resolved component indicator. 

With this type of orifice plate, an error in the perturbed through-flow will result if 
periodic flexture of the plate, analogous to oil canning, takes place. In order to mini- 
mize this problem, a gridwork support structure was attached to the back of the plate. 

Based on the preceding conditions an orifice plate was fabricated (fig. 3) which has 
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1004 holes of 1. 70x10 -meter (0. 067-in. ) diameter and which with a pressure drop of 

O 

345x10 newton per square meter (50 psi) will pass a mean flow of 0. 0631 cubic meter 
per second (2. 23 ft /sec). The multihole orifice plate is held between a pair of NASA 
Orifice Plate Flanges which are located approximately 1. 2 meters (4 ft) from the exit of 


the p um p volute (figs. 1 and 4). The perturbed pressures are measured with two flush- 
diaphragm pressure sensors mounted on the bottom of the outlet line, one on each side 
of the orifice plate (fig. 5). 


Isolation Tanks 

The function of the inlet and outlet isolation tanks is to prevent perturbations of the 
fluid flow from circumnavigating the test loop and to act as system pressure and fluid 
level control systems. These tanks are each approximately 2. 4 meters (8 ft) high and 
1. 2 meters (4 ft) in diameter. Each tank contains a nitrogen-gas-filled rubber bladder 
which fills the upper half of the tank. 


Inlet Flow and Pressure Perturbing System 

The function of this system is to generate the pressure and flow pulsations entering 
the test pump. The system is located between the inlet isolation tank and the beginning 
of the straight inlet line, as shown in figures 1 and 6. The system is composed of two 
parts: the perturbation throttle and the bypass loop. The function of the perturbation 
throttle is to generate the periodic fluctuations to the mean flow approaching the test 
pump. The function of the bypass loop is to allow sufficient flow to bypass the perturba- 
tion throttle to maintain the mean flow required by the test pump operating condition. 

The electro-hydraulic servo-controlled perturbation throttle (fig. 7), which was de- 
signed and fabricated at NASA Lewis Research Center, consists of two multislotted 
plates, one stationary and the other moving. The flow area, and hence the flow, is 
varied by changes in the alinement of the slots in the moving and stationary plates. The 
moving plate is suspended upstream of the stationary plate by four flexure arms which 
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maintain a separation of 5. 08x10 meter (0. 002 in. ) between the plates, thereby elim- 
inating sliding friction while limiting flow leakage. The moving plate is attached to the 
hydraulic drive piston shaft which actuates it back and forth across the stationary plate 
in response to the electro-hydraulic servovalve output. A schematic of the throttle con- 
trol loop is shown in figure 8. 

The throttle as designed has the following characteristics: 

-3 

(1) A full operating stroke of 4. 8x10 meter (0. 19 in. ), which results in low accel- 

eration forces and low stroke velocity 

(2) The ability to operate at full stroke over a frequency range of 1 to 60 hertz 

(3) A total moving weight of 3. 3 newtons (0. 75 lb), resulting in low acceleration 

forces 
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(4) A full open flow area of 18. 5x10 square meter (2. 86 in. ) and an approach 

area of 77. 4xl0~ 4 square meter (12 in. 2 ), giving a high ratio of approach area 
to throat area which results in a linear flow relation with respect to position 

(5) Approximately balanced, lateral, water pressure forces on the moving plate and 

its suspension arms, which result in near symmetrical loading on the hydraulic 
drive piston mechanism 

(6) An actuating piston area of 1. 3xl0~ 4 square meter (0. 2 in. 2 ), which with the 

O 

hydraulic supply pressure of 13 789xlCr newtons per square meter (2000 psi) 
results in a ratio of stall force to moving mass weight of 534 g's (The method 
of valve actuation results in a dc hysteresis of less than 0. 5 percent or 
2. 45X10" 5 meter (0. 001 in. ). ) 

Several other methods of producing the periodic flow fluctuations were considered 
but each was found lacking in one or more characteristics. Other perturbing devices 
considered included the piston perturber, the siren valve, and a hydraulically operated 
commercial valve. The possibility of physically vibrating the pump was also considered. 


Inlet Line Vibration Absorber 

Experience in the line dynamics work (ref. 4) indicates that, when using the acous- 
tical wave equation to measure the perturbed flow exiting from a long line, the longitu- 
dinal motion of the line must be taken into account. To check the accuracy of the meas- 
urements made of the longitudinal vibration velocity of the test facility inlet line, a 
Frahm Vibration Absorber was installed on the line. Discussions on its principle of 
operation are given in references 8 and 9. 

The vibration absorber installed on the facility, shown in figure 9, consists of a 
mass and spring assembly. The mass is a 3. 6-meter ( 13 -ft) long section of 0. 254- 
meter (10-in.) diameter pipe '’weighing” approximately 295 kilograms (650 lbf). The 
mass pipe is split in half along its length and mounted coaxial with the inlet pipe on 
linear ball bearings, and is attached to the inlet pipe by eight specially designed spring 
elements. The linear ball bearing and a typical spring element are shown in detail in 
figure 9. Because of the requirement that the absorber be effective over a large fre- 
quency range, two sets of spring elements of the same configuration were used, where 
the spring constant of one type of element was 1^ times the spring constant of the other. 
For proper operation of the absorber, the spring force must be coaxial with the mass. 
To satisfy this requirement, pairs of diametrically opposite spring elements of like 
spring constant must be connected. By proper combination of spring elements from 
both sets, eight specific frequencies covering the range of 23 to 52 hertz could be set. 
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Instrumentation 


Two types of instrumentation were used: one to make dynamic measurements, and 
another to make steady-state measurements. 

Dynamic measurements were made of the pressure perturbations in the test loop 
and of the longitudinal motion of the inlet line. The dynamic pressure measurements 
were made with commercial, flush-diaphragm, strain-gage pressure sensors. Ten 
dynamic pressures were measured: five in the inlet line at stations 1 to 5; and three at 
the pump, one at the pump inlet and two in the pump scroll (stations 6, 7, and 8, re- 
spectively). The other two dynamic pressure measurements were made at either side 
of the outlet orifice plate (stations 9 and 10). Locations of pressure sensors are indi- 
cated in figures 1, 5, and 10. Calibration of each sensor with its associated amplifier 
and power supply was made by using static calibrating procedures. A schematic of the 
instrumentation signal path and control loop is shown in figure 8. 

The longitudinal motion of the inlet line was measured by using two types of motion 
sensors: a vibration transducer type and a magnet-core - wire-coil type, both units 
located as shown in figure 10. Both units operate on the principle of electromagnetic 
induction due to a magnet moving in a coil of wire. However, their construction and 
method of attachment to the moving item differ. The vibration-transducer-type sensor 
uses an inertial mass as the reference location, and therefore requires a calibration of 
the sensor for phase angle over the test frequency range. The magnet-core - wire-coil- 
type sensor has the coil mounted on a fixed platform and the core mounted on the flange. 
This setup gives output signal voltages proportional to the velocity of the moving core. 
For this sensor the manufacturer’s calibration was verified; however, it is applicable 
only when the core is accurately located in the center of the coil. Because of necessary 
facility maintenance, this condition could not be maintained; therefore, calibration fac- 
tors were determined by using the vibration transducer sensor each time the facility was 
disassembled. 

The output from the pressure and motion sensors is in the form of an alternating- 
current electrical signal. The alternating-current signals are amplified and fed into a 
transfer function analyzer which processes and indicates the data, at any desired refer- 
ence frequency, in resolved component form. The analyzer effectively rejects all fre- 
quencies, noise, and distortion other than the desired reference frequency. The ana- 
lyzer consists of two units: (1) an oscillator which supplies the signal for the 
perturbation throttle controller and a four-phase reference voltage for the resolved 
component indicator, and (2) the resolved component indicator which indicates the in- 
phase and quadrature components of the data signal with respect to the reference voltage. 

Steady- state measurements were made of the steady- state pressure and flow in the 
pump and test loop and of the rotational speed of the test impeller. The pressure meas- 
urements were made with cavity-type, strain-gage, bridge, differential pressure 
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sensors; each pressure sensor was calibrated by using its associated power supply and 
a shared switched amplifier. Measurements were made of the pump impeller inlet total 
pressure, the pump impeller outlet total pressure, and the scroll wall static pressure 
at the locations indicated in figure 10. Both inlet and outlet total pressure were meas- 
ured using a claw-type, total-pressure probe set at the mid-passage location. The inlet 
total-pressure probe was set to the inlet flow angle. The impeller outlet total-pressure 
probe is mounted in an automatic probe actuator and is alined into the flow stream at 
each operating point. Both of these pressure measurements were made relative to 
atmosphere. The wall static pressure in the scroll was measured as a differential with 
the inlet wall static pressure. 

The pressure drop across the outlet orifice plate was measured with a cavity-type 
pressure sensor and with a precision Bourdon-tube-type gage. Both of these measure- 
ments employ pressure taps built into the NASA Orifice Plate Flanges (locations shown 
in fig. 5). 

The mean flow rate was measured with a venturi flowmeter located as shown in fig- 
ure 1. The pressure drop across the venturi meter was measured with an inverted bi- 
fluid manometer using MS 5606 hydraulic oil over water. Both the venturi meter and 
the multihole orifice plate were calibrated in the system using a single-hole, sharp-edge 
orifice plate which was located midway along the length of the 18. 3-meter (60-ft) inlet 
line. The single-hole orifice plate was designed to ASME flowmeter standards. 

The rotational speed of the pump was measured with a magnetic proximity gage and 
a 60-tooth steel gear attached to the pump drive shaft. The indication of pump speed 
was made on an electronic counter. 

METHOD OF ANALYSIS 

Meaningful results will be obtained from this test facility only if reliable values of 
the dynamic flows entering and leaving the pump can be determined. This section dis- 
cusses the methods that are used to make and check these dynamic flow measurements. 
Details on much of the material in this section are presented in reference 4. For the 
analysis and equations used in this experimental investigation, it is assumed that ana- 
lytical stations 1 to 5 are fixed with respect to ground and do not move with the inlet pipe 
and that stations 6, 6A, 9, and 10 are attached to the pipe (see fig. 1 for station loca- 
tions). 
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Inlet Flow Calculation 


To obtain the flow at station 5 of the inlet line with respect to ground, the classical 
one-dimensional wave equation is used in the following form: 


P„ = P(- cos 
n d 


c ol 


n-5 


+ 1 


^ (q 5 ) sin 
Sc V 'n 


c ol 


n-5 


(4) 


Using pressures measured at stations 1 to 4 and at station 5, four independent calcula- 
tions of the flow at station 5 can be made, from which an average Qg flow is 

calculated. n 

At a frequency where sin col /a is near zero for a particular calculation, large 
errors of flow may result from small errors in the measured quantities. At these fre- 
quencies the average flow was determined by omitting that particular calculation. 

To calculate the flow entering the pump with respect to the pump inlet, the velocity 
of the pump inlet and pipe combination must be subtracted from the velocity of the fluid. 
The flow-continuity relation is 


^6A 

S 5 



(5) 


or 


OeA ” Q 5 V 6A S 5 

where Qg^ is the flow leaving the 0. 20-meter (8-in.) pipe or entering the pump, Qg is 
the flow measured with respect to a fixed reference system, and Vg^Sg is the equiva- 
lent vibration flow. 

A check of the measurement of the equivalent vibration flow can be made by deter- 
mining the mechanical admittance of the inlet line. To determine the mechanical ad- 
mittance from experimental data, the driving force tending to vibrate the line in the 
longitudinal direction has to be obtained. In determining the driving force the following 
assumptions are made: 

(1) Momentum forces resulting from flow direction changes are neglected. 

(2) The pressure force at the downstream end of the inlet line is PgSg. 

(3) The pressure force at the upstream end of the inlet line is P-^Sg, and pressure 
forces upstream of the end of the 18. 3 -meter (60-ft) straight pipe are neglected or 
assumed to be counterbalanced. 
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The driving force in the flow direction from station 1 to 5 is then 


F = (P 5 - Pj)®5 


( 6 ) 


For an ideal viscous -damped, single-lumped spring-mass system, the mechanical im- 
pedance is 


Z m =I=D + ] (U-S) 


( 7 ) 


Hence, the mechanical admittance, which is the reciprocal of the mechanical impedance, 


is 


Y = Y = 
m F 


x 


K 


(Xi 


D + j wM - - D - j wM 


D - j Lm - ^ 

*) 


Y m = 


D - j (ujM - — 

(X> 


D 2 + (ooM - -- 

00 


( 8 ) 


which yields at resonance oo = <x Q = yK/M 


Y = Y =I 
m F D 


(8a) 


The mass and spring constant can be determined as follows: 


at resonance 


dZ 


m 


da) 


M + 

2 

oo 


AZ. 


m 


Aw 




oo. 


(8b) 


(8c) 
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Hence, knowing AZ m /Aa> | at co = u> Q , the mass and spring constant of the system 
can be determined from- the experimental data. The mass and spring constant of the 
system can also be determined from static tests. Based on the preceding, two indica- 
tors of the quality of the measurements of pressure and velocity (equivalent vibration 
flow) are 

(1) The agreement of the experimental dynamic values of mechanical admittance of 

the inlet line with the theoretical values determined by using equation (8), 
where the constants D, M, and K are determined by using equations (7) 
and (8c) 

(2) The agreement of the experimental values of mass and spring constant with the 

statically determined values of the same parameters 


Outlet Flow Calculation 

The outlet flow is determined by measuring the pressure difference across the mul- 
tihole orifice plate. This pressure difference is composed of two parts: the pressure 
drop associated with the energy loss due mainly to the sudden expansion at the orifice 
exit, and the pressure difference required to accelerate the mass of fluid between the 
two measuring stations. The pressure drop due to energy loss is proportional to the 
flow rate through the orifice and hence is proportional with respect to the orifice, if a 
linearized equation is used; that is, R Qr is the slope of the mean pressure drop against 
mean flow rate. The pressure difference due to acceleration is proportional to the flow 
with respect to a fixed reference system. Hence, the pressure difference from loca- 
tion 9 to location 10 can be expressed as 


P 9 " P 10 “ Q or R or + J wL 9-10l 


H>r + V 9 S 9^ 


Q, 


Q, 


or 


or 


(9) 


where 


L 9-10 


pi 


9-10 


is the expression for the inertance between locations 9 and 10. 

If VgSg is small with respect to Q or , the solution for Q Qr is simple. Fortu- 
nately, the mechanical vibration of the outlet pipe is small because of the following: 
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(1) The facility structure is considered to be stiffer in the outlet direction than in 
the inlet direction. 

(2) Due to the short length of outlet pipe, the pipe mass is low. This low mass in 
combination with item (1) gives rise to a high resonant frequency and hence to high me- 
chanical impedance at the operating frequencies. 

(3) Because the outlet line is short, the contained fluid mass is low; and therefore 
the driving accelerative pressure forces are low. The low pressure force and high 
mechanical impedance result in a low vibration velocity. Therefore, the outlet flow was 
calculated by using equation (9), where Vg was assumed to be zero, resulting in the 
following equation: 

P 9 “ P 10 = ^or^ R or + i Ct,L, 9_lo^ ( 9a ) 

To decrease the effect of random data scatter in calculating Q or , the pressure dif- 
ference Pg - Pjg was obtained by averaging two sets of data. For one set Pg and 
P^q were measured independently, and for the other set the differential pressure was 
measured directly. 

A check of the quality of the outlet flow measurement is indicated by the agreement 
of the experimentally determined and heuristic values of outlet line inertance, L g_ tank 
and L 10 _ tank . 

The experimental values of outlet line inertance are obtained from the following 
equations: 


9, tot 


Q, 


= u>(L 


9-10 + L 10-tank ) 


or 


10, tot 

Q 

^or 


= wL 


10-tank 


( 10 ) 


(ID 


where it is assumed 

(1) That pipe vibration is zero 

(2) That the numerators on the left side of the equations are exactly Pg ^ - P tank 

and P^g tot - P tank , respectively, but for this analysis P tank is equal to zero 

(3) That Pg , ot and P 10 tot are d y namic total pressures 

Dynamic total pressure is defined as the sum of the measured wall static pressure 
and the kinetic pressure, and is expressed by the following relation which neglects 
second-order terms: 
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(12) 


P tot " P S + 


pQQ 


The heuristic values of outlet line inertance are obtained from the following equa- 
tions: 



where the second term in the equation is the end correction. This end correction ac- 
counts for the small reactance at the end of the outlet pipe "looking" into the outlet tank. 
The end correction is that for a piston mounted in an infinite wall radiating into free 
space and is called the Rayleigh end correction. The end correction is an approxima- 
tion because the physics of the internal assembly of the outlet tank introduces complexi- 
ties which prevent obtaining a more exact expression. A discussion of the limitations to 
the derivation of exact expressions is given in the text of reference 10. The end correc- 
tion is nominally 10 percent of the total inertance. 


Compliance Flow 

The difference between the measured inlet and outlet flow, expressed as 

Q diff = % A ' Q or (15) 

can be caused either by an error in measurement or by a compliance in the pump. 

If the difference in flows is due to a compliance, the following equation, which 
neglects inertance, gives the relation of the compliance, the compliance flow, and the 
perturbed total pressure: 


P 6, tot _ 1 

Q c jwc 


(16) 
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In this equation it is assumed that a single lumped compliance is acted on by a pressure 
located nearest the pump inlet. It is recognized that in the actual case the volume 
change is due to distributed compliance and pressure, but for simplicity the aforemen- 
tioned conditions will be used. The compliance is the ratio of a volume change to an in- 
crease in pressure, where the change in volume is composed of the expansion of the 
pump case and the compression of the fluid associated with the pump. The pressure at 
the pump inlet was chosen because it is the obvious choice if the pump is cavitating and 
a reasonable choice if the pump is not cavitating. The total pressure was arbitrarily 
chosen instead of the static pressure. The difference between the two was less than 
5 percent, or 2° in phase angle. 

A general distinction between measuring errors and compliance flow can be made by 
substituting from equation (15) for Q c in equation (16). The degree of agreement 

of the experimental values of over the frequency range with the theoretical calcu- 

lations of phase angle and magnitude of Q c indicates the degree to which the measured 
difference between inlet and outlet flow is due to compliance rather than to inaccuracies 
in flow measurement. The amount of compliance is determined from the mean value ob- 
tained from equation (16) using for Q c > A check of the value of this compliance 

can be obtained from a static test, in which the volume change in the pump for a given 
increase in pressure is measured. 


TEST RESULTS AND DISCUSSION 

The previous section discussed methods used to indicate the reliability of the per- 
turbed flow values as obtained from measurements of perturbed pressure and inlet pipe 
movement. This section presents and interprets the measured data and calculated pa- 
rameters that are used 

(1) To establish the general reliability of the dynamic pressures, flows, and pipe 

movement 

(2) To demonstrate operation of certain components 

(3) To indicate compliance within the pump 


Test Conditions 

The tests to check the performance of this test facility were made under the follow- 
ing conditions: The frequency of the perturbations to the pump inlet flow and pressure 
covered a range of 2 to 60 hertz. The water in the system was deaerated to an air con- 
tent level of 3 ppm by weight and was maintained at 25. 55° C (78° F). The test pump 
impeller was operated at a rotational speed of 50 revolutions per second. A suction 
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pressure of 482. 63x10 newtons per square meter (70 psig) was maintained while the 
impeller was operating over the mean flow rate range of 0. 0299 to 0. 0524 cubic meters 

Q 

per second (1. 058 to 1. 849 ft /sec). At these conditions, the range of the perturbed 

-4 -4 

flow approaching the pump inlet was 2. 26x10 to 12. 17x10 cubic meters per second 
(0.008 to 0.043 ft 3 /sec). 


Pressure Perturbation Generation and Deduction 

The function of the perturbation throttle was to produce a sinusoidal variation of 
pressure in the inlet line. The relation of the throttle position and the oscillator driving 
voltage with respect to time is shown in figure 11, which is a composite of oscilloscope 
trace photographs covering the entire frequency range over which the throttle was nor- 
mally operated. In each photograph the sweep rate of both traces is the same; this al- 
lows phase comparison of the driving voltage and valve position traces. The oscillator 
driving voltage input to the throttle controller and the sensitivity factors for each trace 
are constant for the series of photographs. The throttle stroke is representative of the 
stroke used during the tests. 

The traces indicate that the throttle is capable of following a sinusoidal driving sig- 
nal with slight distortion. They also show that as the frequency is increased the stroke 
amplitude and phase lag increase. For example, at 60 hertz the stroke amplitude is 
31 percent greater than that at 1 hertz, and the phase lag between the throttle position 
and oscillator driving voltage is 112°, or approximately 1. 8° per hertz. 

A composite of oscilloscope photographs indicating the signals from the flush- 
diaphragm pressure sensors is shown in figure 12. The number by each trace indicates 
the station location of the sensor in the test facility as shown in figure 1. The photo- 
graphs were taken when the pump was in a highly cavitating steady- state condition, 
where the inlet pressure was 68. 947x10 newtons per square meter (10 psig). The fre- 
quency of the pressure fluctuation is 20 hertz, and the amplitude scale is the same for 
all traces. The source of the synchronization signal for all these traces is the oscilla- 
tor signal; therefore, it is possible to compare the phase relation of the traces. For 
example, the pressure at station 5 lags the pressure at station 1 by approximately 90°. 

All the traces, especially those for stations 1 to 4, indicate a general sinusoidal 
shape as a result of throttle movement and system response. 

In the traces of stations 5 and 6 the peaked upper portions and shallow lower por- 
tions indicate the bottoming effect of the cavitation at the pump inlet. This cavitation 
exhibits the characteristic of increasing compliance with decreasing inlet pressure, 
where the effect of the cavitation sets the lower limit on the value of the pressures at 
stations 5 and 6. If the cavitation did not have this characteristic, the shape of traces 
for stations 5 and 6 would be more sinusoidal. Also the traces for stations 1 to 4 would 


16 



be cleaner because of the absence of the reflected waves from the cavitation. This 
cleaning up of the traces was observed in the noncavitating runs. 


Inlet Line Vibration 

A photograph of the oscilloscope traces for the vibration transducer and magnetic- 
core - wire-coil combination is shown in figure 13. The steady-state operating point, 
the dynamic conditions, and the time base for these traces are the same as for those 
shown in figure 12. It is indicated in figure 13 that the line vibration approximates a 
sinusoid of the same frequency as the pressure fluctuations generated in the inlet line. 
Analysis of the traces for the two instruments shows a phase angle difference of approx- 
imately 20° at 20 hertz. This difference is a phase error which is an inherent charac- 
teristic of the vibration transducer and agrees with the value obtained in the calibration 
of the sensor. The accuracy of the vibration velocity as measured by these sensors is 
discussed in a following section. 


Measurement of Perturbed Flows 

Inlet flow . - The quality of the perturbed flow measurement at station 5 in the inlet 
line can be assessed from the plots shown in figures 14 and 15. The ratio of each of the 
four perturbed flows that can be computed at station 5 to the average of these flows at 
station 5 is plotted in figure 14. The figure indicates that the bulk of the flow data lies 
within ±3 percent of the average flow value over the test frequency range. The devia- 
tion of the individual phase angles from the average flow phasor phase angle for the 
same data is plotted in figure 15. This figure indicates that the bulk of the flows has a 
deviation band of ±2° from the average flow angle. It can be inferred from these data 
that the average flow (as defined in the section METHOD OF ANALYSIS) at any frequency 
will have a tolerance equal to or less than ±2° in angle and +3 percent in flow. In com- 
puting the average values for flow and phase angle, the individual data points marked 
with tails were not included; for these points, sin coZ/a is near zero and they are omit- 
ted for the reason discussed in the section METHOD OF ANALYSIS. 

V ibration flow (mechanical admittance) . - The ability to determine the flow due to 
the longitudinal vibration of the inlet line can be deduced from the data in figure 16. The 
figure is a plot of the mechanical admittance of the inlet line, both with and without the 
vibration absorber attached, as a function of the perturbation frequency. The mechan- 
ical admittance is the ratio of the longitudinal velocity to the driving force. The plot 
shows good agreement between the values of mechanical admittance determined experi- 
mentally and those determined using a theoretical relation (eq. (8)). The mass and 
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spring constant used in equation (8) were obtained from experimental dynamic data and 
equation (8c). 

An additional indication of the reliability of the vibration flow measurements is pro- 
vided by a comparison of the values for the mass and the spring constant of the test fa- 
cility inlet structure as determined by two methods: (1) using the experimental dynamic 
data and equation (8c), and (2) performing static pull tests and weight determinations of 
the facility components. 

The mass values were determined to be 920. 875 kilograms (63. 1 slugs) for the ex- 
perimental dynamic data and 980. 710 kilograms (67. 2 slugs) for the weight determina- 
tion. The weight determination mass value is an average of two calculations of the 
mass: one included all the facility components from the inlet tank to and including the 
pump, and the other included only the components from station 1 to and including the 
pump. This was done because there is an uncertainty as to which of the components 
from the tank to station 1 vibrate with the same phase and amplitude as the components 
from station 1 to the pump. 

5 

The spring constant values were determined to be 420. 288x10 newtons per meter 

k 

(240 000 lb/in. ) for the experimental dynamic data and 520. 00x10 newtons per meter 
(297 000 lb/in. ) ±12 percent for the static pull test. The static -pull-test spring constant 
value is an average of several tests in which data were taken in a complete loading and 
unloading cycle. 

The resonant frequency calculated by using the static-pull-test spring constant and 
the calculated weight determinations gives a value of 36. 7 hertz, which is in good agree- 
ment with the indicated resonant frequency of 34 hertz shown in figure 16. The experi- 
mental data at this frequency indicate an average longitudinal movement of the pipe of 
only 7.62x10”^ meter (0.003 in.) rms. 

The agreement of the plots of experimental and theoretical values of mechanical ad- 
mittance and also the agreement between the mass and spring constant values as deter- 
mined by two methods indicate the validity of the vibration measurements. Since the 
equivalent vibration flow is equal to V 6A S 5> it follows that the calculated vibration flow 
is also valid. 

Data shown in figure 16 also illustrate the effectiveness of the vibration absorber to 
reduce movement of the inlet line. The dashed lines through the solid data points indi- 
cate that these points were taken with a specific combination of springs whose design 
resonant frequency is indicated on the line. The data with the absorber attached were 
taken at the same steady-state and perturbed operating conditions as the data of run 1 
without the absorber. The marked reduction of admittance with the absorber attached 
is an indication that the longitudinal vibration of the line was substantially reduced. The 
frequency span over which a particular combination of spring elements was used is in- 
dicated by the dashed curves through the data points. The rapid rise from a minimum 
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value in a small frequency span indicates that each spring element combination is 
sharply timed. 

Outlet flow . - The perturbed flow leaving the pump is determined by using equa- 
tion (9a), which assumes that the vibration velocity of the outlet line is zero and is 
therefore neglected. An indication of the validity of the outlet perturbed flow measure- 
ment can be obtained by comparing the experimental values of outlet line inertance (eqs. 
(10) and (11)) with the heuristic values (eqs. (13) and (14)). This comparison is given in 
figure 17 for several noncavitating test runs. Data below 12 hertz have not been in- 
cluded because of the large data scatter caused by the extremely small pressure at 
station 10. These plots indicate that the experimental values fall into a band which ap- 
proximates the heuristic values up to 28 hertz. In this frequency interval the mean ex- 
perimental value is 3 percent higher than the heutristic value. Above 28 hertz the data 
band rises to a higher value and remains at this value to 52 hertz. The shift in level of 
the experimental data above 28 hertz is probably due to a change in the radiation impe- 
dance at the pipe-to-tank junction. See section METHOD OF ANALYSIS for a discussion 
of this radiation impedance. In both of the plots shown, the limits of the data band are 
about ±5 percent of the mean experimental value and remain the same over the entire 
frequency range. 

The consistency of the experimental values for both measuring stations and the good 
agreement of the mean values to the respective heuristic values indicates that the meas- 
urements of perturbed outlet flow are reliable. 

Combined flows . - The relation of the inlet, outlet, and vibration flows is illustrated 
in figures 18 and 19. The magnitude of these flows over the range of perturbation fre- 
quency is shown in figure 18. Because these flows are phasor quantities, they cannot be 
manipulated arithmetically. In figure 18 the first point of interest is the concurrence of 
the flows into and out of the pump over the entire frequency range. For almost the en- 
tire range there is a good match, indicating that the method of accounting for the vibra- 
tion flow works well. The second point of interest is the magnitude and trend of the 
vibration flow, which at low frequencies is very small and has little effect on the flow 
determination into the pump. At a frequency of 18 hertz the vibration flow starts to in- 
crease with increasing frequency until its magnitude is equal to the flow calculated at 
station 5 in the line. Thus, the vibration flow has a marked effect on the determination 
of the flow entering the pump. It is obvious that if there is no compensation for this vi- 
bration flow at frequencies above 18 hertz, determination of the flows into the test pump 
will be in error. The third point of interest is that over the entire frequency spectrum 
the outlet flow is consistently larger than the inlet flow. This is significant in that it 
indicates the presence of a form of storage or compliance within the pump. 

The magnitude and phase relation between the various flows is shown in figure 19. 
The distance from the origin to any point is the magnitude of the flow measured. The 
angle measured clockwise from the right horizontal abscissa to the magnitude line is the 
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lag angle of the flow relative to the master oscillator signal. The flow phasor which is 
the line between the station 5 flow and the pump inlet flow point is the vibration flow. 

The flow phasor which is the line between the outlet flow point and inlet flow point is the 
difference flow. 

Further corroboration of the validity of the inlet line flow, equivalent vibration flow, 
and pump outlet flow is presented in figures 19 to 21. It is indicated in figure 19 that the 
difference flow Q^yy is obtained by the phasor subtraction of the equivalent vibration 
flow and the pump outlet flow from the flow at station 5. Because of the requirement of 
phasor subtraction, any large error in any of these three flows will show up as an in- 
consistency in the resultant value of Q^y*. The error which does exist in the values of 
these three flow measurements can be evaluated with figures 20 and 21. A plot of the 
reactive portion of the difference flow impedance with perturbation frequency is shown 
in figure 20. The data are from seven test rims covering the steady- state noncavitating 
operating range of the test pump. The plot shows that the data follow the general trend 
of a straight line (labeled theoretical curve) with a slope of -1. The plot also indicates 
that the majority of the data points fall within a band of ±33 percent about the value of 
the theoretical straight line. The significance of this data band is indicated by the plot 
of the ratio of the difference flow to the summation of the magnitude of the three meas- 
ured flows with perturbation frequency shown in figure 21. The data in this figure are 
from the same test runs as used for figure 20. The plot of figure 21 indicates that the 
data have systematic variation over the frequency range, with a nominal value of 6 per- 
cent. Based on this mean value of 6 percent, a worst case of ±2 percent change in the 
measured flows will result in a ±33 percent change in the value of Q^yy (i. e. , 

0.02/0.06) and, similarly, in a ±33 percent change in the value of jPg/Q^ycl ’ 
is the value of the data band shown in figure 20. ' ^ 

The variation of the data with frequency in figure 21 is interesting in that, when the 
flow ratio is at its lowest values, the largest deviation appears in figure 20. Also shown 
in figure 20 in the 25- to 35-hertz region is an interesting deviation from the general 
trend of the data; the deviation is indicated by the dashed line. The exact reason for its 
existence is not known, but it is speculated that the fluctuation is associated with a res- 
onance of the inlet line and the front face structure of the pump. 

Compliance flow . - The phase angle between the difference flow and outlet flow for 
the data in figure 19 suggests that the difference flow is a compliance flow and not a 
gross error in measurements. This is further substantiated by the information contained 
in figures 20 and 22. The slope of -1 for the straight line labeled theoretical curve in 
figure 20 exhibits the theoretical slope of a compliance as defined by equation (16) when 
plotted on log-log paper. This suggests that Q^.y in equation (15) is a compliance flow. 
This is corroborated by the data in figure 22, which is a plot of the phase angle of the 
difference flow impedance with perturbed frequency. The plot shows that the data gen- 
erally lie within ±10° of -90°, which is the general characteristic of a compliance. 
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Since has both of the general characteristics of a compliance, it can be concluded 

that is largely a compliance flow. 

In order to make an independent evaluation of the existence and magnitude of the 

compliance, static -pressure tests were performed on the test pump. The objective of 

the tests was to determine the compliance of both the pump structure and the volume of 

water it contained. The data obtained in these static tests give a value for the total 

compliance of 182. 802x10"^ meter"* per newton (3. 091x10 ftVlb). This value is for 

the pump structure and only the water contained in the pump. If the compliance of the 

volume of water from the exit of the pump to the orifice plate is included, the value of 

the total theoretical compliance is 300. 490x10" ^ meter ® per newton (5. 081x10"^ ft^/lh). 

Both of these values are close to and bracket the value of the pump compliance deter- 

“13 5 “8 

mined from perturbed data, which is 265. 953x10 meter per newton (4. 497x10 
ft^/lb). This suggests that the difference between the perturbed pump inlet and outlet 
flows is largely a compliance flow. A difference between the static and perturbed test 
values of compliance is to be expected because the exact perturbed test conditions could 
not be duplicated in the static tests; that is, the pressure distributions in the pump and 
the volumes of water were not the same for the static and perturbed tests. Included in 
the plots of figures 20 and 22 are data points (solid circles) from a test run which had 
the vibration absorber attached to the test facility. The good agreement of the values of 
compliance impedance and its phase angle for the data with the vibration absorber at- 
tached confirms conclusions previously arrived at from the following facts: 

(1) The vibration flow measurements under nonabsorbed operation are valid. 

(2) The difference flow is actually a compliance flow and not an error in measure- 
ment. 


CONCLUDING REMARKS 

This report has discussed the general design concepts and preliminary operation of 
a facility designed to make experimental evaluations of pump performance with imposed 
perturbations to the pressure and flow at the pump inlet. The results from the prelim- 
inary operation of the test facility indicate that the necessary measurements required to 
obtain the dynamic performance of the pump can be made. 

The most significant point of the test is that the methods chosen to measure the 
three basic perturbed flows yield valid results. The methods involve the measurement 
of the perturbed pressure drop across a special multihole orifice plate to measure the 
flow leaving the pump, the use of the acoustical wave equation in a long hydraulic line to 
measure the flow approaching the pump, and the measurement of the longitudinal veloc- 
ity of the inlet line to determine the equivalent vibrational flow. The difference of the 
latter two flows gives the flow entering the pump. 
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It is indicated that for noncavitating pump operation the difference between the per- 
turbed flow into and out of the pump is a compliance flow. The value of the largest ex- 

_7 

perimental noncavitating compliance determined, a volume change of 6. 56x10 cubic 
meters (0. 04 in. ), was corroborated by the data obtained in static tests of the compli- 
ance of the liquid in the pump and the dimensional change in the pump structure. 

The data indicate that the longitudinal motion of the inlet line is very important and 
must be taken into account when using the acoustical wave equation to measure perturbed 
flows. This was found to be especially important in the region of the point of mechanical 
resonance of the inlet system. In this region the vibration flow can equal the calculated 
line flow and hence produce large errors in the evaluation of the dynamic characteristics 
of the pump if it is not accounted for. 

The tests show that the specially designed perturbation throttle was capable of fol- 
lowing a sinusoidal driving signal and producing a usable sinusoidal variation of pressure 
over the frequency range of 2 to 60 hertz. 

The preliminary results gained by using the Frahm Vibration Absorber indicate 
that this type of device can be used to minimize the longitudinal vibration of the test fa- 
cility. However, this type of device has a narrow frequency range for effective opera- 
tion. 

Based on the experience gained in these tests, careful consideration should be given 
to the placement of the dynamic pressure sensors in the pump. The pressure sensors 
should be placed where flow areas can be specifically defined; this will aid in the calcu- 
lation of the velocity head portion of the total perturbed pressure. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 13, 1971, 

720-03. 
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APPENDIX A 


SYMBOLS 


A pump gain 

a wave velocity, m/sec (ft/sec) 

B pump impedance 

c compliance, m^/N (ft^/lb) 

D effective mechanical viscous-damping resistance, N-sec/m (lb-sec/ft) 

F vibromotive force, N (lb) rms, complex 

f frequency, Hz 

o 

j defined by j = - 1 

K effective mechanical stiffness, N/m (lb/ft) 

L inertance, N-sec^/m^ (lb-sec^)/ft^ 

l length, m (ft) 

M mass, kg (slug) 

o o 

P pressure perturbation, N/m (lb/ft ) rms, complex 

o o 

Q volume flow perturbation, m /sec (ft /sec) rms, complex 

— 3 3 

Q volume flow steady state, m /sec (ft /sec) 

3 3 

Qdiff defined by eq. (15), m /sec (ft /sec) rms, complex 

5 5 

R resistance, N-sec/m (Ib-sec/ft ) 

2 2 

S cross-sectional flow area, m (ft ) 

V vibration velocity, m/sec (ft/sec) rms, complex 

Y m mechanical admittance, m/N-sec (ft/lb- sec), complex 

c c 

Z impedance, N-sec/m ° (lb-sec/ft 0 ) 

Z m mechanical impedance, F/V, N-sec/m (lb-sec/ft), complex 
p density, kg/m^ (slugs/ft^) 

a) 2ni 

Subscripts: 

c compliance (flow), see eq. (16) 

d discharge 
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n location number (1 to 4 as applicable) 
o resonance 

or orifice 

s inlet 

tank outlet tank 

tot total (calculated from the sum of the measured static and the computed dynamic), 
see eq. (12) 

vib vibration 

X real part of a complex number 
Y imaginary part of a complex number 

1-10 locations shown in fig. 1 (Analytical stations 1-5 are attached to ground; analyt- 
ical stations 6-10 are attached to the pump and pipe. ) 

6A see fig. 1 
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APPENDIX B 


EXPERIMENTAL DETERMINATION OF ACOUSTIC VELOCITY 

The acoustic velocity of a liquid in a round uniform pipe of known dimensions is 
easily computed. Because the long inlet pipe in this facility was made up of rolled 
welded pipe which has normal dimensional variations, an experimental determination of 
the acoustic velocity of the liquid in this line was made. The acoustic velocity was de- 
termined by using a version of the classical Kundt's method. This method involves the 
solution of the acoustical wave equation (1) to yield the frequency corresponding to a 
half wavelength in the line. At 


^ 1-5 

a 


27tU 1-5 
a 


equation (1) reduces to 


P 



(Bl) 


Putting this in complex notation gives 


Equating real and imaginary parts, 


P 1Y " P 5X ‘ J P 5Y 

(B2) 

IX = " P 5X 

(B3) 

1Y = " P 5Y 

(B4) 


Multiplying equation (B3) by equation (B4) yields 


P 1X P 5Y " P 1Y P 5X 


Putting these pressure values in terms of their indicator values yields 

g l T X E lX g 5 T Y E 5Y = g l T Y E lY g 5 T X E 5X 


(B5) 


(B6) 
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where 


g pickup sensitivity 
T indicator sensitivity 

E indicator reading 

Dividing equation (B6) by g^^g^Ty yields 

E 1X E 5Y = E 1Y E 5X 


(B7) 


The products E^E^y anc * E ^yE obtained from several runs were plotted as 

a function of frequency. The frequency at which the difference of the products goes to 
zero is used to calculate the acoustic velocity of the liquid in the line. The average 
value of the experimental acoustic velocity was found to be 1195 meters per second 
(3919 ft/sec) with an average difference between runs of 5. 8 meters per second 
(19 ft/sec), or 1/2 percent. The theoretical values of acoustic velocity for the dimen- 
sional tolerance range of the rolled steel pipe yield values of 1143 to 1207 meters per 
second (3749 to 3961 ft/sec), which encompass the experimentally determined value. 

The simplicity and accuracy of this experimental method, which does not require 
calibration of either the sensors or indicating apparatus, makes this a desirable means 
of determining the acoustic velocity of a liquid in a pipeline. 
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Figure 1. - Pump Perturbations Test Facility, showing analytical reference stations. For analytical purposes stations lto 5 and 6A are 
assumed to be attached to ground; stations 6 to 10 are assumed to be attached to the pump and pipe. 
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Figure 10. - Test pump and instrumentation. 





















(c) Frequency, 40 hertz. (d) Frequency, 60 hertz. 


Figure 11. - Comparison of oscillator drive voltage and valve position. Valve stroke (peak to peak) equals 18.75 percent of full stroke 
(peak to peak) at 1 hertz. Top trace indicates valve position; bottom trace indicates drive voltage. 
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Ratio of individually calculated flows at station 5 to average of these flows 









Difference between phase angle of individual flows at station 5 
and average phase angle, deg 






Mechanical admittance of inlet line, 1/Z_, m/N-sec 


Orifice impedance 



Figure 16. - Mechanical admittance of inlet line as function of perturbation frequency. Pump inlet pressure for 
all runs, 482. 63x 10 3 newtons per square meter (70 psi). 
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Figure 17. - Outlet line inertance as function of perturbation frequency. Pump inlet pressure for all runs, 
482. 63xl0 3 newtons per square meter (70 psil. 
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Figure 18. - Magnitude of perturbed flow as function of perturbation frequency. 
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Figure 19. - Typical example of perturbed flow in polar coordinates. 
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Figure 20. - Reactance portion of difference flow impedance as function of per- 
turbation frequency. Pump inlet pressure for all runs, 482. 63x 10 3 newtons per 
square meter (70 psi). 
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Figure 21. - Ratio of difference flow to summation of measured flows as function of perturbation 
frequency. Pump inlet pressure for all runs, 482. 63x10^ newtons per square meter (70 psi). 
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Figure 22. - Phase angle of difference flow impedance as function of perturbation frequency. Pump 
inlet pressure for all runs, 482.63x10^ newtons per square meter (70 psi). 
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